Cosmic distance-duality as probe of exotic physics and acceleration 
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In cosmology, distances based on standard candles (e.g. supernovae) and standard rulers (e.g. 
baryon oscillations) agree as long as three conditions are met: (1) photon number is conserved, 
(2) gravity is described by a metric theory with (3) photons travelling on unique null geodesies. 
This is the content of distance-duality (the reciprocity relation) which can be violated by exotic 
physics. Here we analyse the implications of the latest cosmological data sets for distance-duality. 
While broadly in agreement and confirming acceleration we find a 2-sigma violation caused by excess 
brightening of SNIa at z > 0.5, perhaps due to lensing magnification bias. This brightening has 
been interpreted as evidence for a late-time transition in the dark energy but because it is not seen 
in the oIa data we argue against such an interpretation. Our results do, however, rule out significant 
SNIa evolution and extinction: the "replenishing" grey-dust model with no cosmic acceleration is 
excluded at more than 4-sigma despite this being the best-fit to SNIa data alone, thereby illustrating 
the power of distance-duality even with current data sets. 



I. INTRODUCTION 

In 1933 Etherington 

HUH proved a beautiful and 
general duality that implies that distances in cosmology 
based on a metric theory of gravity are unique: whether 
one uses the apparent luminosity of standard candles 
(yielding the luminosity distance, dr,(z)) or the appar- 
ent size of standard rulers (the angular-diameter distance 
(Ia(z)), does not matter since they are linked by distance- 
duality jjtfj: 
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where z is redshift. Distance-duality holds for general 
metric theories of gravity in any background (not just 
FLRW) in which photons travel on unique null geodesies 
and is essentially equivalent to Liouville's theorem in ki- 
netic theory. While it is impervious to gravitational lens- 
ing (for infinitesimal geodesic bundles) it depends cru- 
cially on photon conservation. Our aim in this paper is to 
discuss how distance-duality may become a powerful test 
of a wide range of both exotic and fairly mundane physics 
and to present a general analysis of what constraints on 
violations of distance-duality arise from current data as 
well as critically analysing the conclusions drawn from 
recent type-la supernovae data (also discussed in the 
appendix) . 

To test distance-duality we use the latest type la su- 
pernovae (SNIa) data 0, |E IE as a measure of the 
luminosity distance, dt(z) @|. This data includes a sig- 
nificant number of z > 1 observations. Our estimates of 
the angular-diameter distance, d A {z), come from FRIIb 
radio galaxies 0, lioll. co mpact radio sources 0, 0, Il3| 
and X-ray clusters 14] . It is important to remember 
that some of this data predated the discovery of accel- 
eration by SNIa and that there are now completely in- 
dependent, indirect, estimates of d A , e.g. from analysis 



of the 2QZ quasar survey [H (giving n A = 0.71±g;??) 
and strong lensing from a combination of the CLASS 
and SDSS surveys with a maximum likelihood value of 
Ha = 0.74 — 0.78 m good agreement with estimates 
from radio sources. 

All these data sets broadly agree with an accelerat- 
ing, high-f^A cosmology. Nevertheless, there are a few 
observations in disagreement with the accelerating 'con- 
cordance' model (e.g. |17|'). there are suggestions that 
SNIa may suffer from significant extinction ^8|, evolu- 
tion 01 °r axion-photon mixing |2l|. There are also 
radical alternatives to general relativity, such as MOND 
[20j . Distance-duality gives us a way to test all of these 
possibilities. 



II. DISTANCE-DUALITY VIOLATIONS 

Since our aim in this paper is to promote distance- 
duality as a powerful test of fundamental physics it seems 
appropriate to begin by describing some phenomena that 
could be detected through violations of distance-duality. 
The most radical violations would arise from deviations 
from a metric theory of gravity or in cases where photons 
do not travel on (unique) null geodesies (e.g. torsion 
or birefringence). Other interesting possibilities include 
variation of fundamental constants such as G, but we do 
not discuss any of these possibilities here because they 
are either already tightly constrained, difficult to give 
predictions for or are too implausible given current prior 
beliefs about gravity. 

Instead we restrict ourselves to phenomena that may 
reasonably occur given our understanding of particle 
physics or astrophysics. We will also not discuss obvi- 
ous possible sources of violation such as unaccounted for 
systematic error or biases in estimates of either dj, and 
d A . While this would be obvious first place to look to 
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explain a violation of distance- duality we consider it to 
be trivial and hence will not discuss it further except 
when we put limits on the size of such effects later using 
distance-duality. 



A. Photon number violation 

Perhaps one of the most likely sources of duality- 
violation is non-conservation of photon number. This 
could have a mundane origin (scattering from dust or 
free electrons) or an exotic origin (e.g. photon decay or 
photon mixing with other li ght states such as the dark 
energy, dilaton or axion |U|23|). However, all of these 
effects tend to reduce the number of photons in a light 
bundle and therefore reduce the apparent luminosity of 
a source. If unaccounted for, this dimming makes the 
source appear more distant, i.e. increases cLl- Since (1a 
is typically unaffected (or negligibly altered) by such ef- 
fects, this rather generally implies that the ratio in equa- 
tion (1) becomes greater than unity. The case of axion- 
photon mixing has been studied in [22| and the results 
there show that this type of dimming cannot obviate the 
need for cosmic acceleration. 

We can parametrise scattering or loss of photons by 
studying the photon Boltzmann equation integrated over 
frequency allowing for a collision functional: 



(2) 



n 7 is the number density of photons and 7, a are con- 
stants that control the scattering/decay cross-section of 
the photon. Ho is the current value of the Hubble con- 
stant, a = —2 corresponds to a scattering cross-section 
oc Pcdm oc (1 + z) 3 , as in the case of Compton scatter- 
ing from free-electrons. The case of photon decay corre- 
sponds to a — 1. 7 > implies loss of photons. In fact 
7^0 leads to a violation of distance-duality that grows 
roughly exponentially with rcdshift, see Fig (Q. 

For the sake of generality we also consider power-law 
deformations of distance-duality that parametrise our ig- 
norance about the effects of more exotic physics. This 
yields a 3-parameter (a,f3, 7) extension of equation 
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where E(z) = H(z)/Hq is the dimensionless Hubble ex- 
pansion normalised to unity today. Distance-duality cor- 
responds to (/3, 7) = (1, 0) (in which case a is arbitrary). 



B. Lensing and finite beams 

Distance-duality holds exactly only for infinitesimal 
light bundles in which case gravitational lensing has 
no effect on the duality. In practice however, observa- 
tions are made with different finite-sized bundles. Es- 
timates based on observations on large angular scales, 
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FIG. 1: Graphic evidence for violation of distance- 
duality. The binned data for dL(z) (triangles, SNIa) and 
dA(z) (circles) are shown in equivalent magnitudes relative 
to the flat concordance model (Qa = 0.7, Q m = 0.3) with la 
error bars. They should coincide if distance-duality holds but 
they differ significantly at z > 0.7. The dashed curves are the 
best-fit FLRW models to the dji(z) data (top) and dz,(z) (bot- 
tom) separately with no loss of photons (7 = 0). The solid 
curves have the same underlying FLRW model (Qa = 0.81, 
fim = 0.22) but the lower curve includes the best-fit bright- 
ening (7 = —0.036, see eq. 2) with a = —2,f3 = 1. Since 
the violation of distance-duality increases exponentially when 
7/O more high redshift data and/or smaller error bars will 
significantly improve constraints. 



(such as the 2df 10QZ survey [l5( or the proposed KAOS 
survey 38] ) will be very weakly affected by gravitational 
lensing, while SNIa observations may be strongly affected 
by lensing (by an amount up to 0.3 mag [UISI or more), 
depending on the fraction of compact objects in the uni- 
verse. Using such different techniques to estimate dA 
and c?l implies that lensing will violate distance-duality 
by an amount that depends on the fraction of compact 
objects [2j| . This opens the interesting possibility that 
future data will be able to test the fraction of compact 
objects by searching for such lensing-induced violations 
of distance-duality. 

One way to get around this lensing-induced violation is 
to analyse objects that can give both tii, and dA- An in- 
teresting possibility in this category are type 2 SN where 
dA can be estimated from observations of the photo- 
sphere. Unfortunately dA data of this sort is currently 
limited to very low redshift [25| . 
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FIG. 2: Supernovae are brighter relative to d,A data. 

7-Aa likelihood plot in the case a = —2, /3 = 1 which corre- 
sponds to a photon scattering probability oc (1 + z) 3 . The 
best fit corresponds to 7 = —0.036, i.e. brightening of SNIa 
relative to the d,A data, as required from figure 1. The very ex- 
tended, diagonal, contours are the weak 1 and 2a constraints 
found using only the SNIa data. This illustrates the power 
of blending cLl and cLa data as a consistency check of existing 
data and as a test of new physics. 



III. CONSTRAINTS FROM CURRENT DATA 

Here we use the standard FLRW equations to calculate 
the theoretical distance oIa(z) as a function of the cosmic 
parameters (Omj^a) (over which we then marginalise, 
as they are determined by the angular diameter distance 
data) and use © to infer cIl{z) given (a, (3, 7). 

We use a standard Markov-Chain Monte Carlo method 
with one chain of 10 6 points per model to sample the like- 
lihood and derive the marginalised limits. It is probably 
useful to point out here that the 3-parameter eq. |31 con- 
tains two "artificial" degeneracies clearly visible in the 
likelihood contours of Fig © and Q: firstly if 7 = 
then a is completely unconstrained (and /3 f» 1) and sec- 
ondly there is a value of a around for which the integral 
is very close to a logarithm, and so the full right hand 
side of the equation becomes approximately (1 + 2;) /3+7_1 
which leads to the strong degeneracy visible in fig. 0] 
More details of our data sets and method are given in 
appendix A. 

For this reason it is preferable to study the absolute 
goodness of the fit in the full three-dimensional parame- 
ter space of (a, (3, 7). Hence instead of marginalising we 
show plots found by maximising (equivalent to marginal- 
isation in the case of a Gaussian likelihood). In this way, 
it is easier to see where the well-fitting models are lo- 
cated. When quoting the limits on 7 we do of course 
marginalise. 

In Figure (1) we show the binned cIl(z) and (1a(z) data 
as a function of redshift converted to magnitude (relative 



to the flat concordance model) assuming distance-duality 
holds, in which case both data sets should lie on the same 
curve. The shaded regions shows the effect of the best-fit 
7 = —0.036 (a = —2,13 = 1) on the underlying oIa{z) 
showing how it is possible to simultaneously fit the oIl 
and d,A data with a single model. Also shown are the very 
different best-fits to the cLl and cLa data taken separately. 
While the <1a data favour a flat universe, the SNIa data 
favour a very closed model (ruled out from the CMB) 
due to the unexpected brightening at z > 0.5. 

Figures (3) & (4) show the results of our Monte-Carlo 
Markov Chain likelihood analysis. It is clear that the 
distance-duality prediction {[3 = 1 , 7 = 0) is not favoured 
by current data with the best fit occurring in the de- 
generacy region at the edge of the figure, {a, (3, 7) = 
(0.1,4.0,-2.7) with = 217. In comparison, the 

best-fit FLRW model, (/?, 7 ) = (1,0), has xLin = 223 - 

Figure (2) shows the joint 7 — fl\ likelihood that fol- 
lows when one imposes (3—1 and a — — 2 by assuming 
scattering from objects whose number density scales as 
(1 + z) 3 (such as Compton scattering by free-electrons) 
we find that the best fit for the absorption coefficient 
is 7 = -0.036 with xLin = 219 and -0.07 < 7 < 
at 95% confidence. Surprisingly, the best-fit corresponds 
not to absorption but brightening, as is clear from fig- 
ure (1) since the <1a(z) data lies above the di, points. 
The extra parameters are justified in all cases from the 
Akaikc information criterion while the Bayesian infor- 
mation criterion favours introducing 7 at the expense of 
curvature, while it marginally disfavours introducing all 
of (a, /?, 7) depending on the binning of the data (see 
e.g. |3(j). The magnitude of the effect corresponds to 
an increase of about 5% in the number of photons per 
Hubble time, a very large violation of photon conser- 
vation. We can put this into perspective by compar- 
ing it with the expected loss of photons due to Comp- 
ton scattering by the free electrons in the ionised inter- 
galactic medium. At z < 3 helium is expected to be 
doubly ionised (fy = 0.5), leading to a free-electron den- 
sity n e — QbPcritfi — YfY)/rriM where Y = 0.24 is the 
primordial helium abundance. We therefore find a scat- 
tering amplitude of "(Compton = o- T n e /(2H ) ~ 10~ 3 , a 
factor of about 50 less than the best-fit (and of opposite 
sign). 

A plausible explanation is magnification bias through 
gravitational lensing. If distant SNIa are preferentially 
detected if they are brightened then this would cause an 
apparent violation of the reciprocity relation as discussed 
in section (|II Bll . It has recently been demonstrated that 
lensing does significantly affect current high-z SNIa sam- 
ples [26J. Brighter high-z SNIa are preferentially found 
behind overdense regions of galaxies and can differ from 
demagnified SNIa by 0.3 — OAmag. The induced bias 
may be sufficient to provide the ~ 0.1 — 0.2mag bright- 
ening required to remove the violation of distance-duality 
we have documented above. Alternatively, since smaller 
compact radio sources are typically brighter |l3j| , in an in- 
complete magnitude limited survey, high-redshift sources 
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FIG. 3: Goodness-of-fit contours (1 and 2 a) for the 
parameters a and f3 of equation (J3J. Distance-duality im- 
plies (3—1. If photons decay at a constant rate, then a = 1. 
Higher values of a correspond to models where the photon 
scattering cross section increases as the universe expands. If 
photons are affected by "scatterers" with a constant comoving 
density, then a = —2. a > 1 corresponds to a rather unphysi- 
cal region of parameter space where the probability of photons 
scattering increases with the expansion of the universe. 



will be systematically smaller and yield a larger value of 
d A {z). 



A. Ruling out replenishing dust 

Riess et al Q found that the best-fit model to all cur- 
rently available SNIa was not an accelerating ACDM 
model but rather a replenishing grey-dust model |27| 
with A = which causes redshift-dependent dimming 
of the SNIa, with a changing from —2 to 1 at z = 0.5. If 
this was the correct explanation then we should expect 
a marked violation of distance duality with the (La data 
lying below the di, data since it would correspond to a 
non-accelerating universe. Our results show that this is 
not the case (indeed we have the opposite problem!) 

A detailed analysis of this model based on H?],!^ gives 
a best-fit to all the data off^A = 0.77±0. 13 showing that 
the combined data, in contrast to the SNIa data alone, 
rule out the replenishing dust model at over 4-cr. 



B. Is dark energy evolving? 

While we have discussed interesting physics which 
violates distance-duality, dark energy dynamics is not 
among them. Hence evidence from SNIa for significant 
evolution in the dark energy equation of state w(z) at low 
redshift 0, I2J H3 now appear less significant since 
the signal is not seen in the c?a(z) data. As our mea- 
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FIG. 4: Goodness-of-fit contours (1 and 2 a) for the 
parameters a and (3 of equation Distance-duality 
implies /3 = 1 and 7 = 0, which corresponds to photon con- 
servation. This point is acceptable for the three-parameter 
case, but becomes unfavoured when we limit ourselves to the 
sub-space (/? = 1, a = —2) of constant comoving-density scat- 
terers. 



surements of distance-duality improve we will be able to 
obtain better constraints on dark energy evolution. 



IV. CONCLUSIONS 

In this paper we have emphasised distance-duality as 
a test of fundamental and exotic physics related to the 
metric nature of gravity and photon conservation on cos- 
mic scales. Although stringent constraints will arise in 
the next few years the test is already proving power- 
ful. In particular we are able to essentially rule-out non- 
accelerating models of the universe which explain the su- 
pernova dimming by grey-dust scattering, extinction or 
evolution. Interestingly, current data suggests a small 
(2-er) discrepancy that may be due to lcnsing-induccd 
magnification bias of the high-z SNIa. 

One can ask if it will be possible to distinguish vio- 
lation of photon conservation from deviations of grav- 
ity from a metric theory assuming systematics errors are 
eliminated. One interesting way to do this would be to 
use binary black holes as standard gravitational wave 
candles [34j to give an independent estimate of dr,(z). 
Comparing this against the dr,(z) found from SNIa us- 
ing e.g. the JDEM/SNAP satellite and against distance- 
duality should allow us to distinguish between the two 
possibilities. Further, large galaxy surveys such as the 
proposed KAOS experiment, will provide accurate esti- 
mates of cU(z) out to z — 3 35], allowing us to test devi- 
ations from distance-duality at the level of a few percent, 
implying that this diagnostic will mature into a unique 
and powerful test of fundamental physics on cosmological 
scales. 
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APPENDIX A: DATA SET DETAILS 

The main supernova data set (cLa) is the "gold" subset 
of Riess et al [4| . We checked that this gives essentially 
the same results as the earlier data in Tonry et al. j5|, 
Barris et al. @, and Knop et al. QUI (TBK). Although 
we used our own code to evaluate the resulting likelihood, 
it follows closely the one of John Tonry and gives the 
same results. 

For the estimates we used the data sets of Dal y fc 
Djorgovski (DD), Gurvits [HI (G) and Jackson |3 
(J). DD provide their data directly as dimensionless y(z) 
and we use it in this form. G gives the data as angular 
sizes 9{z) with oIa — 1/9 and we need to marginalise over 
the unknown "standard ruler" I. This is analogous to the 
case of supernovae. For this reason the radio galaxy data 



also does not depend on the Hubble constant. J also pro- 
vides angular sizes, but pre-binned and uses error bars 
determined so that the resulting \ 2 value per degree of 
freedom is unity. We then marginalise over an indepen- 
dent angular size in this case as well. We checked that 
we obtain the same confidence regions as [T^ | when using 
the J data set alone. 

How stable are our results to changes of the underlying 
data sets? Taking the absorption model as a test case, 
leaving out any single d,A data set does not change the 
constraints on 7 appreciably, although if we drop DD 
then 7 = becomes acceptable at 2a. If on the other 
hand we use the combined supernova data sets of TBK, 
we find stronger evidence for a violation of distance- 
duality, with 7 < —0.01 at 2a. The conclusion that there 
is something systematically different between the SNIa 
and the RG data sets is therefore rather stable. 

As a further test of the radio galaxy data, we have 
included the gas mass fraction data of Allen et al., and 
marginalised over all nuisance parameters (in this case 
the bias, the Hubble constant and the baryon density). 
The X-ray data is consistent with the RG data, and its 
addition does not change the constraints on 7. But as 
eq. (15) in |14| is only given for flat universes, we quote 
our results without this data set. Strong constraints may 
also come in the future from SZ data l36l. 
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